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detren, in a static, unidirectional, inhomogeneous magnetic field with finite 



' flux can be calculated from the massive Euclidean Schwinger model's determi- 

nant, detach) in the same field by integrating detach, over the fermion's mass. 



Since det^en for general fields is central to QED, it is desirable to have nonper- 



' turbative information on this determinant, even for the restricted magnetic 

' fields considered here. To this end we continue our study of the physically 

relevant determinant detach- It is shown that the contribution of the massless 
' Schwinger model to det sch is cancelled by a contribution from the massive 



sector of QED in 1 + 1 dimensions and that zero modes are suppressed in 

det Sch- We then calculate det sch analytically in the presence of a finite flux, 

cylindrical magnetic field. Its behaviour for large flux and small fermion 

^ • mass suggests that the zero-energy bound states of the two-dimensional Pauli 

, 

Hamiltonian are the controlling factor in the growth of In detach- Evidence 
is presented that det sch does not converge to the determinant of the mass- 
less Schwinger model in the small mass limit for finite, nonzero fiux magnetic 
fields. 
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I. INTRODUCTION 



Fermion determinants produce an effective measure for the boson fields of a Euclidean 
field theory when the fermions are integrated out, as shown long ago by Matthews and Salam 
[1]. Such determinants in an external boson field, or a random boson field with a cutoff, 
are infinite dimensional and need to be defined [2] before the boson fields can be integrated. 
Once defined, their analysis is notoriously difficult, especially if they possess a symmetry 
that should be preserved. For example, most classical estimates of fermion determinants [3] 
invariant under a local U(l) transformation violate this invariance. Lack of nonperturbative 
information on fermion determinants is reflected in the necessity to make loop expansions or 
the more extreme quenched (valence) approximation. As a result physical effects predicted 
by the theory may be lost. 

Nonperturbative information on the fermion determinant, such as its growth in the com- 
plex coupling plane, is central to an analysis of the nature of the perturbation series of the 
associated field theory [4] . Intuition tells one that Fermi statistics, visible in the alternating 
signs of the determinant's loop expansion, ought to slow down the growth of a perturbation 
series with order. 

There is also the question of stability. Specifically, in the case of an Abelian gauge field, 
the measure is Gaussian, so that if the fermion determinant grows faster than an inverted 
Gaussian, it is doubtful that it is integrable respect to the gauge field's measure. 

Given the sparseness of nonperturbative information on physically relevant fermion de- 
terminants in four dimensions we thought it useful to try to find a solvable example for a 
broad class of boson fields. An obvious choice is the fermion determinant of quantum elec- 
trodynamics in time-independent, unidirectional, inhomogeneous magnetic fields. Although 
of physical interest, this choice of fields suffers from the fact that they are a set of measure 
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zero as far as the functional integral over the vector potential is concerned. Nevertheless, the 
fermion determinant for such magnetic fields remains unsolved except for the special case of 
a homogeneous magnetic field that was dealt with over half a century ago by Euler, Heisen- 
berg, and Weisskopf [5] and later on again by Schwinger [6] . A thorough understanding of 
this problem would be helpful for a more general understanding of the physical content of 
fermion determinants in quantum electrodynamics. As we will see below, there are some 
significant simplifications that recommend this problem for analysis. 

This paper is organised as follows. In Sec. II we review previous relevant results. Section 
III justifies our restricted choice of magnetic fields on which the results of Sec. II rely. We 
then go on in Sec. IV to discuss the suppression of zero modes in the fermion determinant 
of the massive Schwinger model. As Sec. II makes clear, the four-dimensional determinant 
is obtained by integrating this determinant over the fermion's mass. Section IV also illus- 
trates the profound change in the two-dimensional determinant when the fermions are given 
a mass. In Sec. V we calculate the massive Schwinger model's determinant for a finite, 
nonzero fiux magnetic field. Section VI contains a discussion of the zero-mass limit of the 
massive Schwinger model's determinant, and Sec. VII summarizes our results. 



In a previous paper [7] it was shown that the renormalized fermion determinant, det^en, 
of quantum electrodynamics in four dimensions (QED4) in a smooth, static, unidirectional 
magnetic field with fast decrease at infinity is related to the Euclidean, massive Schwinger 
model's (Euclidean QED2's) determinant, detsch, in the same magnetic field. Specifically, 



II. REVIEW OF PREVIOUS RESULTS 
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where the gauge invariant dets is related to detsch by 

Since the first term in (2.2) is the contribution to the effective action from the second- 
order vacuum polarization graph, In dets may be viewed as the sum of all one-loop fermion 
graphs in two dimensions beginning in fourth order. Our designation of this sum by In dets 
follows the definition of Seller [2] and Simon [3], with the third-order graph vanishing by 
C-invariance. Equation (2.1) states that the sum of the corresponding graphs in QED4 is 
obtained by integrating In dets over the fermion mass m. The first term in (2.1) is the contri- 
bution to the effective action of QED4 from the on-shell renormalized second-order vacuum 
polarization graph. The function B is the Fourier transform of the magnetic field, which 
may be asumed to point along the 2;-axis, in which case V\\ is the volume of the zt box. In 
both (2.1) and (2.2) the determinants are defined by Schwinger's proper time definition [6]. 
Note that the charge, e, will always occur in the combination eB in position space, which 
has the invariant dimension L~^. Note also that in Euclidean QED2 potentials associated 
with unidirectional magnetic fields are a set of measure one. 

The lesson of (2.1) is that the massive Schwinger model [8] is more than a model in view 
of its direct bearing on physics in four dimensions. Unlike the determinant of the massless 
Schwinger model [9], det ^c/i in (2-2) is not known explicitly. Nevertheless there are some 
important results. One of these is an expression of the paramagnetism of fermions in an 
external magnetic field as summarized by the "diamagnetic bound" 

detsch < 1, (2.3) 

for > 0[2, 10, 11]. This bound also expresses the positivity of the one-loop effective 
action for Euclidean QED2. As a consequence of (2.3) we were able to put an upper bound 
on the growth of det^en in QED4 for strong, unidirectional, inhomogeneous magnetic fields. 
Performing the dilation — > XA^ on the vector potential and letting A become large we 
obtained 
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where = y" d\B\r) [7]. 

There are other results for detsch imphcit in the hterature. For example, if falls off 
sufficiently rapidly so that j (fr\A^\'^ < oo for all g > |, then one can relate detschie) to 
the zeros of the determinant, considered as a function of a complex coupling e [12]. It is 
known [13] that for m 7^ these zeros occur in quartets e„, — e„, e* and — e* , and therefore 
detsch cannot vanish for real e for these potentials [14]. 



We will defer the discussion of the result of Haba [15] for det^en until Sec.III. 



In view of the direct connection between QED4 and QED2 in the case of unidirectional 
magnetic fields it was thought worthwhile to obtain more specific information on det sch 
that would enable one to make use of (2.1). In a second paper [16] it was shown that the 
exact calculation of detsch reduces to a problem in nonrelativistic, supersymmetric quantum 
mechanics. That is to say, 

Indetsc/. = ^ / (iV09V + dXTr{[(H+ + m')"^ - (H_ + m')-», (2.5) 

ZTT J Jo 

where the supersymmetric operator pair H± = (p — AA)^ ^ XB are obtained from the two- 
dimensional Pauli Hamiltonian in (3.3) below. The auxiliary potential, (p, is related to the 
vector potential by A^^ — €^,^8^(1) and to the magnetic field by B — —d^(f) or 

4){r) = -— / dVln |r - r'\B{r'). (2.6) 
27r J 

The antisymmetric tensor e^u is normalized to €±2 — 1. Again, the starting point for the 
derivation of (2.5) was the proper time definition of det5c?i, which respects gauge invariance 
and allows one to select the Lorentz gauge. Our representation of detsch makes a sharp 



5 



separation between the contribution from the massless Schwinger model, the first term in 
(2.5), and its massive counterpart. We have not integrated the first term in (2.5) by parts 
as is usually done. For nonzero flux fields, falls off like ^, and therefore an integration 
by parts is invalid in this case. 



Within the Lorentz gauge there is still the restricted gauge freedom of — > + c, where 
c is a constant. Since detsch is gauge invariant, the term proportional to c in (2.5) must 
vanish, in which case 



27r 



2m^ / dXTr 



{H+ + mY'-iH^ + my\, (2.7) 

since —9^0 = B and $ = j cfrB. Differentiating (2.7) with respect to e gives the index 
theorem on a two-dimensional Euclidean manifold [17], 



— = rn^Tr 
2n 



(H+ + mY-(H_ + m^) 
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n+-n_ + -^[5^0)-5L(0)], 
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(2.8) 



where n± denote the number of positive and negative chirality zero-energy bound states of 
H±, and (^±(0) are the zero-energy phase shifts for scattering by the Hamiltonians H± in a 
suitable angular momentum basis I. Thus, the index theorem in QED2 follows from gauge 
invariance. 



Now suppose (2.5) is written as 

Indetsch^^ [ d'^r(f)d'^(f> + 2 r dXTr{[Pi^\Q) - P^^\Q)](f)} + nonzero modes, (2.9) 
2tt J Jo 

where Pi^\o) are the zero-mode spectral measures associated with H±. Noting that 

— i7_ ~ B and assuming B sufficiently weak, we calculated the second term in (2.9) in 
first Born approximation and showed that it cancelled the massless Schwinger model term. 
We conjectured that this was true in general, and we will show that it is in Sec. IV. 
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Finally, we wish to retract a claim in [16]. By the Aharonov-Casher theorem [18], n+(n_) 
e|$| 

, depending on whether e$ > 0(e$ < 0), where [x] stands for the nearest 



are given by 



27r 



|e$| 

integer less than x and [0] = 0. Thus, if > 1 there are zero-energy bound states, and it 

2ti 

was stated that when these are included in the second term of (2.9) the logarithmic growth 
of 4> and the slow, algebraic fall off of the bound-state wave functions for large r would cause 
detsch to vanish. This is false. In fact, zero modes are suppressed in detsch as we will see in 
Sec. IV. 



III. CHOICE OF FIELDS 

The fermion determinant is part of a functional integral whose measure, dii{A), is that of 
the gauge-fixed, free Maxwell field. As diJ,{A) may be reahzed on <S', the space of tempered 
distributions, we seem to be stuck with rough potentials that are hard to analyze. A way out 
is to realize that there is a logarithmic ultraviolet divergence in QED2 due to the vacuum 
energy graph shown in Fig.l that has to be regularized and subtracted out. One way to 
regularize [19] is to smooth by convoluting it with an ultraviolet cutoff function h^^S, the 
functions of rapid decrease. That is, let A^ — A/^* h\ so that A^ is a polynomial bounded 
C°° function, meaning that for each a — (m, n) there is a N{a) and a C{a) with 

<C(l + x^ + yY. (3.1) 

The regularized free photon propagator is 

Jdn{A)A^^{x)At{y) = D^^^{x-y), (3.2) 

whose Fourier transform, D^j^{q), is proportional to |/iA(g)p, where is the Fourier trans- 
form of /iA- One possibility is to choose Ha = with hx{q) = 1 for < and hA{q) = 
for q'^ > 2h?. Note that the measure diJi{A) is not regularized. Thus, without loss of gen- 
erality we may assume that our potentials are smooth and polynomial bounded. Hereafter 
we will drop the superscript A and denote these potentials simply as A^. 



Qm Qr. 



At{x,y) 
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Now in order that the fermion determinant exist it seems necessary to require the mag- 
netic fields derived from the potentials to have finite fiux. The reason for this restriction 
is connected with the degeneracy of the ground state associated with the two-dimensional 
Pauli Hamiltonian 

H ={p-Af- a^B > 0. (3.3) 

Specifically, Avron and Seller [20] have considered the class of polynomial, infinite-fiux mag- 
netic fields 

N 

B{r)^J2Ur-Cnf''\ (3.4) 

n=0 

where {A„} and {C„} are arbitrary real numbers and {kn} and N are nonnegative integers. 
They have shown that the ground state of H is infinitely degenerate and that the mani- 
fold of zero-energy bound state wave functions is parametrized by a point in R,2(2fcmaa;+i)^ 
irrespective of the translational invariance of the magnetic field. In the constant field case, 
N — 0, one has point spectrum, and the vector space is a plane whose points specify the 
center of rotation of Landau orbits, corresponding to the known degeneracy \eB\LxLy/2TT, 
with Lx,Ly — > oo. This degeneracy persists for all excited states in the constant field case, 
and we suspect that the excited states for fields with N > are at least as degenerate as 
j^2(2A;max+i) Although we will show in Sec. IV that zero- modes are suppressed in our repre- 
sentation (2.5) of detsch, it cannot make sense out of the volume-like divergences associated 
with the degeneracy of the excited modes. Nor are we aware of any definition of determinant 
that can. Accordingly, we will confine our attention here to finite flux flelds. In our view, 
a determinant of an infinite fiux field should be considered as the limit, if it exists, of a 
determinant calculated from a field confined to local planar regions. 

Restriction to finite fiux is also consistent with the additional need to regulate volume 
divergences before taking the thermodynamic limit of the Euclidean Green's functions. Such 
divergences appear in the vacuum energy graphs, including the one in Fig.l, when the de- 
terminant is integrated with respect to dii{A). A volume cutoff may be introduced via the 
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determinant, assumed previously calculated for a generic class of fields, by replacing B{r) 
with g{r)B{T), where 5^ is a suitable volume cutoff such as geC^ or g = xa, the characteristic 
function of a bounded region A C R^. 



At this point we mention the result of Haba [15], who studied detj-en in the field F^j^^^x) = 
Gptv + F'nui^), where G/^i, is constant, and is sufficiently smooth and bounded. It was 
concluded that with the introduction of a space-time volume cutoff A on the otherwise 
infinite flux field F^i^ that 

Indet^en = ^{G^ - l\G*G\) ln(l + G') 
-G^ I d^xb{G,F'), 

J A 

where b is bounded in G. This saturates our upper bound (2.4) for the special case — 
and G12 — B. Haba then goes on to use this result to obtain evidence suggesting the 
instability (trivality) of QED4. In view of the potential importance of this result and of 
our foregoing remarks concerning infinite flux fields, it would be worthwhile to repeat the 
calculation with F^,^ confined to a finite region A C from the beginning and repeating 
Haba's estimates with A held fixed. 



As a consequence of the need to regularize and the finite flux condition we may confine 
our attention to smooth, polynomial bounded potentials that, in the Lorentz gauge, have 
the asymptotic form 

Mr) = ~-^ + 0(^,), (3.5) 

where $ is the flux associated with B and r >> a, the range of B. Our long-term goal is to 
study detsch under the scaling B XB, A — > 00 as well as the determinant's mass depen- 
dence. This latter point, as we now know, is especially relevant to determining Indet^en in 
QED4 for unidirectional, static, inhomogenous magnetic fields. 
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We know of no previous calculation of a finite flux determinant associated with massive 
fermions in two or more dimensions. In Sec. V we will calculate such a determinant analyt- 
ically for a magnetic field confined to a thin cylindrical shell with radius a: 



ZTT a 

Although B is not derivable from a polynomial bounded potential it has the virtue that the 
(5(r)/r-type singularity of a magnetic fiux string is absent and that detsch exists for this 
field. It is an instructive example, especially as we believe it gives an insight to the matters 
raised above. 



IV. SUPPRESSION ZERO MODES IN detgch 
A. Conventions 

Consider a Dirac fermion in a static, unidirectional magnetic field directed along the 
z-axis. Its Hamiltonian in the xy plane is 



H = 7°7.(p - eA) + 7°m, 



(4.1) 



with the 7 matricesQy" = a^, 7^^ = —ia'', k = 1,2. It has the structure 



H 



m L 



V 



(4.2) 



L'^ —m 

where L is a linear differential operator, and is the Hermitian conjugate of L. The 
positive- and negative-energy eigenfunctions of H, 



^Where possible we adopt the notation and conventions of Jaroszewicz [21] whose analysis of the 
chiral anomaly associated with the Hamiltonian (4.1) in a solenoidal magnetic field is relevant to 
the work presented here. 
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Hi;EAr)^Ei;EAr), (4.3) 

are normalized to 

I dV4,,(r)V'B, ,,(r) = 5^y5(E - E'), (4.4) 

in the energy continuum, where A is a degeneracy parameter. As a consequence of (4.2) the 
eigenvalues satisfy E'^ >-rv?. Since 

= (p - eKf - eBa^ + m^ (4.5) 

we have 

r,r,t 1 

= (p - eAf T eB, (4.6) 

so that the supersymmetric operator pair H± in Sec. Ill are given by — LL\ if_ = L^L. 
The eigenfunctions of and if_ will be denoted by u^ji and v^j^i, respectively, 

H+Uyj^i{v) = wuy,^i{v) 

H_Vyj,i(r) ^wvyj,i(r), (4.7) 

where / is a degeneracy paramater and w — E^ — w?. Their normalization in the energy 
continuum, w > 0, is 

j (frul i{v)u^,^i,{v) = 5{w - w')5ir. (4.8) 

Since the continuum extends down to = ±m or ty = we expect the states on the edge of 
the continuum to show up as unbounded resonances [22]. When |e$|/27r > 1, bound states 
will also appear at the bottom of the continuum [18], and care must be taken to include 
these in the completeness relations for ^l)E,\ and Uw^i^Vu,,i- 
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B. Suppression of Zero Modes 



Referring back to (2.5) consider 

d\udeisch/de = - / dVt/)^^^ 
vr J 



+2m^ J d^r(j){r) {r\{H+ + m^)-^ - {H_ + m^) 
Now consider the large mass limit of 



2\-l 



(4.9) 



m 







poo „ p 

dte-'^'i^BM+Oit)) 

Jo ZTl 



(4.10) 



where we used the heat kernel asymptotic expansion [16] 

l±teB{r) + Q{f) 









1 int L 



(4.11) 



Combining (4.7) and (4.10) we get 



lim ( r 



m 



m 



hm / dw- 



roo 

/ dwY.{\u^,{v)\'-K,{v)\') 

JO I 



roo 

= / {r\dP+(w)-dPJ 
Jo 



w]\r) = 



eB{r) 



(4.12) 



where P+{w) — P-{w) is the difference of spectral measures associated with H± and where 

roo 

H± = / wdP±(w). 
Jo 



Equation (4.12) has the following physical interpretation. The charge density induced in 
the vacuum by the background magnetic field is 



(/(r)) 



^^E(l^i^,A(r)p-|^-i.,A(r)n, 

A 



(4.13) 
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which is determined by the spectral asymmetry density 

r](E,r) = ^discEtr {r\(E - H)-^ - (£; + //)-» 

ZTT 

-Y.i\^EAr)\'-\i^-E,x{r)\'), (4.14) 

A 

where discEh{E) = h{E + ie) — h{E — ie). Due to the structure of H in (4.2), rj may also be 
expressed in terms of the eigenfunctions of H± in (4.7) [21,23] as 

r){E,r) = 2mY,i\uE^-mM{r)\^ - \vE^-mM{r)\'), (4.15) 
I 

so that 

(/(r)) = / dw , E(l^wWI' - lvKr)r). (4.16) 
Comparing (4.16) with (4.12) we see that the limit in (4.12) is equivalent to 

hm {f{r)) = _^!^. (4.17) 

The natural length scale here is the range of the magnetic field, and so (4.17) states that the 
vacuum change density induced by a magnetic field whose range is large compared to the 
fermion's Compton wave length is — e^S(r)/47r, in agreement with a remark by Jaroszewicz 
[21]. 

Now consider (4.9) again in the form 

d\ndetsch/de = - [ d^rcpd'^cp 

+2 / dV0(r) r dw-^ y:{\u^,{v)\' - \v^,{r)\') 

J •/ Jib "x (JJ ^ 

n J 

/poo on 
(fr4>(r) (1 -){r\dPJw)-dPJw)\r). (4.18) 
Jo w + 



Prom the last line of (4.12) we finally obtain 
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d\udeisch/de (fr(l)B + 2 J dV0(eS/27r) 



-2 / d^c 



w 



-2Tr 



w + ni? 
H. H. 



{r\dP+{w) - dP-{w)\r) 



(4.19) 



which shows exphcitly that zero modes are suppressed in detsch- It also shows how substan- 
tially mass has altered the fermionic determinant of QED2: the contribution of the massless 
Schwinger model to detsch has been cancelled by a contribution from the massive sector of 
QED2. 



C. Gauge Invariance of (4.19) 

If we revert back to the definitions (2.5) or (4.9) of detsch, then the invariance of detsch 
under the restricted gauge transformation (f)+c results in the index theorem (2.8). Recall 
that the massless Schwinger model's contribution to detsch was critical to establishing this 
result. Now, in (4.19), we see that the massless Schwinger model's contribution has been 
cancelled by a contribution from the massive sector of det sch that contains the zero modes 
of H±. Invariance of (4.19) under 4> ^ (p + c now requires that 

/ d'r r dw-^Y.iM^)\" - - 0. (4.20) 

J JO UJ "Y" fit ^ 

In the remainder of this section we will recall the known result that the integrand in (4.20) 
can be expressed as the divergence of a current [21,23] and then go on to shown explicitly 
that (4.20) is true in the case of radial symmetry. 

Differentiating (4.7) with respect to w, denoted here by an overdot, one easily obtains 

[21] 

K,i{r)f-K,ii^)\^^V-S^,i{r), for w > 0, (4.21) 

where 
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S^,; = -[K,i^^w,i - C^ul^i)uyj,i - 2iul iAuy,^i -(u^ v)]. (4.22) 

To make further progress we will specialize to the case of radial symmetry so that the degen- 
eracy paramater I is identified with angular momentum. We would then like to interchange 
the space integral with the sum over partial waves in (4.20) in order to convert the space 
integral to an integral over a circle at infinity. Jaroszewicz has already discussed this inter- 
change in another context [21], and we repeat his reasoning here. Since if+ — if_ ~ B, and 
B has finite range, a, the difference between the wave functions u^j^i and v^j^i, for fixed w, 
decreases with increasing I due to the rising centrifugal barrier that excludes them from the 
region where B{r) ^ 0. The energy required for the wave function to penetrate the region 
r < a is of the order of w > ^, suggesting that the partial wave sum in (4.20) is effectively 
cut off at \l\ ~ ^/wa. Accepting this reasoning, we get 

where dl = TiRd9 in the r, 9 plane. Since the potential A may be assumed to be a pure 
gauge field at infinity tangential to 5"^, the A-dependent terms in (4.22) may be dropped. 
Indeed our potentials (3.5) approach vortex fields that manifestly satisfy this assumption. 
Using the asymptotic form of u^^i, 



Tr ( --^^ - ^" = r dw—^Yl lim / S^KR')-^^1, (4-23) 



Uw,i{f-id) ~ 2 ■^t: '^r '^e' 
r a 



cos 



(4.24) 

and similarly for v^^i, where Sf{w) and Sf{w) are the scattering phase shifts, one gets [21] 

lim / iCRYdl =-y \SY(w) - Sr(w)] , for w > 0, (4.25) 

where the overdot continues to denote differentiation with respect to w. The factor e~''^ 
instead of e*'^ in (4.24) is for later notational convenience. Due to the supersymmetry of the 
operator pair H± we have from (4.7), 

L^H+u^^i{r) = H_L''u^^iir) = wL^u^,i{r), (4.26) 

which indicates that L^Uw^i oc Vw^i and hence that 
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Sf{w) = S^_^{w), modTTfor^/; > 0, (4.27) 

in agreement with Jaroszewicz [21]. In deriving (4.27) we used 

= e'\^ -ieAr + -^ + eAe) 
or r od 

■fi , d i d . / . 

r a or r oO 27rr 

where A = {Ar,Ae). The condition for the invariance of (4.19) under the restricted gauge 
transformation ^ + c is now reduced to 

OO L 



L^oo 



= hm [d\(w) - SI, Aw)] = 0, for w > 0. (4.29) 

This is physically reasonable since the wave equations for u^^i and Vw,i become scale invariant 
outside the range of B, where the potentials V± defined by 

H± = -d^ + V±, (4.30) 

have a ^ behaviour and to where u^^i and Vw,i are mainly confined due to the enormous 
centrifugal barrier building up as |/| — > oo. Hence, we do indeed expect lim^^oo — 0- 



V. detsch IN A CYLINDRICAL MAGNETIC FIELD 



A. The Green's Functions 



In this section we will use the representation (4.9) of detsch, as it is simpler than (2.5) 
and more instructive than (4.19). For the magnetic field in (3.6) we have for the auxiliary 
potential given by (2.6), 



(f)(r) = -— / dVln |r - r'\B{r') 
2tt J 



— < 



* 1 

m r, r > a 

2n 

ma, r < a. 

27r 



(5.1) 
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Since detsch is invariant under ^ + c we may let c — — In a and use the potential 

ZTT 



(l>{r) = -—\n^-je{r-a), (5.2) 
in which case the contribution of the massless Schwinger model is eliminated straightaway: 

= 0. (5.3) 

The associated vector potential is 

A = —e(r-a)e 

27rr ^ ^ 

^ ^e. (5.4) 

r 

Referring to (4.9) we see that the Green's functions G±^i defined by 

l=—oo 
-1 oo 

= ^ E G^,{k;ry)e-^^^'-''\ (5.5) 

are central to the calculation of detsch- The radial wave functions of H+^i and i?-,/, denoted 
by Uk2^i{r) and ffe2^/(r), respectively, satisfy 



(i^ 1 (/ + e$(r))2 e$\r; v / x ,2 / x 



f Id (/ + e$(r))2 e$'(r)\ ,x ,9 /x /r 

These equations have hnearly independent solutions H^{kr) for r < a and H^^{kr) for 
r > a, where if J" and H~ denote the Hankel functions Hj^^^ and Hj^^\ respectively. The 
dimensionless constant e$/27r has been denoted simply by $. 

The calculation is simplified by introducing the Green's functions 

g±,i{k;r,r') = V^G±,i{k;r,r'), (5.7) 
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where 



(5.8) 



and 



2 _ 1 

4 



$'(r) 



(5.9) 



The outgoing- wave Green's functions ^^/(C = ± = ±1) are constructed from [24] 



gli{k;r, r') 



(5.10) 



where 4> and /+ are regular and irregular solutions, respectively, of 



(5.11) 



i7+ for ^ = ± are the associated Jost functions, and r<,r> denote the lesser and larger values 
of r, r', respectively. The indices I and ^ have been omitted on the right-hand side of (5.10) 
to reduce notational clutter. 

Regular solutions of (5.11) are 



k l'ly^J/(A;r 



r < a, 
, r> a, 



(5.12) 



[aH+^^{kr) + (3Hf^^{kr) 
and the irregular, outgoing wave solutions are 

^/kr [AH^(kr) + BHf(kr)] , r<a 
y/kr H^_^{kr), r > a. 

The constants ct, /3, A, B are determined by the joining conditions at r = a obtained from 
(5.11) and (5.9): 



(5.13) 



(f){k, a—) = (j){k^ a+) 
(l)'{k, a+) - (l)\k, a-) = -C*0(^, a+)/a, 



(5.14) 



and similarly for /+. These give 



a — -Ti^ka Ji(ka) Hi_^^_^{ka) — Ji^^{ka)Hi^^{ka) 
13 = --K^ka YJi-s,{ka)H^j^i^{ka) - Ji{ka)H^^_^{ka) 
A = -ir^ka H^_i^{ka) Hf_^{ka) — Hf (ka)H,t^.^_^{kc 



B = -n^ka \H^{ka)H^^^_^{ka) - H^_^{ka)H+^{ka) 



(5.15) 



Finally, the Jost functions are 



k--2-V\l3rW[H+^{kr), H[_^^{kr)] 



pk-^ 

TT 



(5.16) 



where W is the Wronskian, and (3 is given in (5.15). 



At this point g^^i is fully determined by (5.10), (5.12), (5.13), (5.15) and (5.16). We 
could now continue ^ into the upper half of the k plane by letting k — me*^, ^ — > f with 
m > 0, thereby making contact with the Green's functions {t\{H± + m^)~^|r') in (4.9). We 
will do this later. But first we want to demonstrate that when the flux is sufficiently large 
for fixed Z, Q^^i does indeed contain a zero-energy bound state. This will be done by deriving 
a completeness relation. 



B. Completeness 

The proof of the completeness of the bound and scattering wave functions associated 
with 'H±^i will follow the general procedure outlined by Newton [24], taking into account 
that contains a ^-type potential. 



Consider the integral 

Jkk dr'h{r')gli{k;r,r') 
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= - / dkk r dr'hir') 
Jc Jo 

/oo 
dr'h{r') 



,^(p{k, r')f+{k,r) 
n(t>{k,r)f+{k,r') 



IC Jr ' ' J+{k) 

= Ii + h, (5.17) 

where h{r) is square-integrable, and C is the contour in Fig. 2. The contribution of /ir to Ii 
from the large semicircle F is evaluated by using the asymptotic forms of 0, /+ and ^7+ in 
the upper k plane. From (5.15) and (5.16), 

J+ = ^aA;^~l'l [ji-^{ka)H,^^{ka) - Ji{ka)H,^^_^{ka) 

^ -'^k"^-\i\e-^-^/\ (5.18) 

From (5.12), (5.13) and (5.15) one gets, in the upper k plane, 

(t){k,r) ~ cos(A;r - - W) 

|A;| — >• oo V TT 

|A;| ^ oo V TT 

and hence for \R\ — oo, 

Iir ^-ij^dk dr'h{r')e'^^'-^''^--^''^ cos{kr' - \ttI - \ti) 
~ - "^hir) j^dk f^dr'e'^^'-''^ 

- -mh{r). (5.20) 

The contribution of I2T to I2 from the contour F may be defined by replacing the upper 
limit of the r' integration in (5.17) by r + and letting // — > 00 later [24]. Then for 

\R\ — > 00, 

~ -i j^dk j"^^^ dr'h{r')e'^^''-'^''^--^''^ cos{kr - ^ttI - ^tt) 
~ -^/i(r) [ dk f^'^dr'e'''^''-'^ 

2 Jv Jr 

~ \mh{r). (5.21) 
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Due to the analytic properties of H^{ka) and Ji{ka) in the upper half of the k plane, J^+{k) 
has no zeros for Imk > 0. For example, on the positive imaginary /c-axis we have from 
(5.18), 

J+{e-^ka) = [Ii_^{ka) Ki+^{ka) + Ii{ka)Ki+^.^{ka)] , (5.22) 

TT 

where Ky and Ii arc modified Bcssel functions. Since Ky[x), Ii{x) are positive for x > 
0, j7+(e*'^/^A;a) is manifestly free of zeros for k > 0. Therefore, 

J^dkkJ dr'h{r')gli{k]r,r') = 0. (5.23) 

Combining (5.17) with (5.20), (5.21) and (5.23) we get 

+ - r^^dr'h(r') [ dkkgti(k;r,r'). (5.24) 

TT Jo J-y 



Equation (5.24) may be simplified by the following relations, vahd for positive k, 

(f){ke''',r)^(f){k,r) 
f+{ke'\r) = -e^"/2e-^'^('+*)/;(A;,r) 
J+ike'"") = e^'^e-'^'^^^e-'^'^J^k), (5.25) 

where 

J_{k)^W{f^,cl>) 

Aia , 1 III , , 

= fc2-KI (5.26) 

71 

f4k, r) = ^/k^Hf^^ikr), r>a, (5.27) 

and where a is given by (5.15). For real k, J^{k) = J-{k). Equations (5.25), (5.26) and the 
result, valid for real k, 

2k 

Im [J4k)U{k, r)] = —(f){k, r), r < a, (5.28) 

TT 

reduce (5.24) to 
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h{r) 



dr 



'h{r') I 



oo 

dk-, 



\Mk)\ 



-(j){k, r)(j){k, r') + - dr'h{r') / dkkgli{k; r, r'). 

TT Jo Jj ' 

(5.29) 



The integral around the semicircle 7 will contribute only if Q^i develops a second-order pole 
at A; = 0. A tedious calculation confirms that this happens for ^ = +, Z < 0, Z + $ > 1 and 
for ^ = —1, / > 0, / + $ < —1. This result is in accord with a remark by Jaroszewicz [21], 
who stated this would happen for $ ^ integer; our result holds for all admissible values of 
The residue of the double pole is the zero-energy bound state ipi of 'H±x. 

lim / dkkgti{k;r;r') = -inipi{r)ijji{r') , (5.30) 



for ^{l + ^) > 1,^1 <Q and where 



2(1 + 


/|)(l 


/4 




-1) 


a2 







a, 



r < a 
r > a. 



(5.31) 



The bound states are normahzed: 



/■oo 

Jo 



1. 



(5.32) 



Combining (5.30) with (5.29) and letting // — > 00 gives the completeness relation 

4 1"°^ k"^ 

i,+ ^^ |j7-^(fe)|2 '^^^' + = - r'). (5.33) 

We have gone through this calculation to show that Q^^^ does indeed contain a bound 
state for sufficiently large flux and how it manifests itself. Later, when we consider the small 
mass limit of det ^c/i, one should keep in mind that this is being controlled by the bound 
states, since the m — > limit is approaching the second-order pole of Q^^i along the positive 
imaginary /c-axis. 



C. Calculation of detgc/i 

The outgoing wave Green's function Q^i{k;r,r') in (5.10) may be continued to A; = 
im, m > 0. According to the definition of detsch in (4.9) we will only need the difference of 
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Q^^i and G-^i, which simphfies the calculation. Due to (5.2), we may confine our attention 
to r = r' > a. Then from (5.10), (5.12), (5.13), (5.16) and (5.26), 



[ J-{k) 




J-ik) 


J 






J+{k) 





gX^i{k]r,r) -g+i{k]r,r) = — 



which is related to the (S'-matrix for the Ith partial wave by 

Mkf ■ 

In (5.34) the Jost ratios can be rearranged to give 

J- _ 2 Ji-({ka)Ji+^{ka) - Ji{ka)Ji+q,_i:{ka) 



W+dkr)f 



- 1. 



J+ Ji_^{ka)Hj^^{ka) - Ji{ka)Hj^^_^{ka) 
We now let k me'^,m> in (5.34), (5.36) and (5.8), using [25] 

J,(ame^'^/') = e^^'^/'7,(ma) 

i/+(ame^"/2) = —e-''"'/^KJma), 



(5.34) 



(5.35) 



(5.36) 



(5.37) 



to obtain 



2\-l 



r) — (r 



2\-l 



— —r 



$-1 



///, 



- 



i^f+$(mr). 



(5.38) 



Jl-lKi+^ + Ii+iKi+q, + IiKi+^^i J 

where all modified Bessel functions in the brackets have argument ma. If the two terms in 
(5.38) are combined by a common denominator then a remarkable simplification occurs to 
give 



= 2$r < [ma 



d 



In 



Lima) 



-1 



dma \Ki_^_<^(ma) 

2 



12 _ $2 



X < 



f Kf_^^[mr 
Kf+^{rna 

Ii{mr) 



-, r > a 



r < a, 



(5.39) 



. I? {ma)'' 

where we have put on record the result for r < a. As a check on our results we verifed that 
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m 



J 7 _„ 



m 



2\-l 



l=—oo 



r ) — (r 



2\-l 



(5.40) 



for > 0, in accordance with (2.8). This was done by interchanging the integral and 

sum, which is allowed since the series in (5.39) is uniformly convergent for r > 0, and then 

defining Y.'i^-oo ^ hniL-^oo Z^^l- The integrals of the modified Bessel functions are known 

[26]. Again, the reader is reminded that $ denotes e$/27r, where $ is the flux, and that this 

symbol is active from (5.6) onward. 

Finally, combining (4.9), (5.2), (5.5), (5.7) and (5.39) gives 

dlndetsch r , ^ ^ Kf+dmar) 
-4 ma $ / drrmr > — — ^ 

d , f lAma) 



X 



2 



ma- 



in 



(5.41) 



dma \Ki^q,(ma) 

The term in the curly brackets is positive for all Figure 3 displays the numerical calcula- 
tion of the right-hand side of (5.41) for the cases ma — 1 and 10~^ for < $ < 999.5. The 
plots were generated for half- integral values of The data in both cases are consistent with 
a logarithmic growth of dlndetsch/d^ with $ given by — ln$-|- const, where the constant 
is about 0.3 for ma — 1 and —2.5 for ma — 10"^, consistent with the "diamagnetic bound" 
in Eq. (2.3). 



The integral in (5.41) can be calculated exphcity for integer Assume this is the case, 
and let ^ = N. Since the sum in (5.41) is uniformly convergent for r > 1, we may let 
I —>■ I — N and interchange sum and integral. Since K_i = Ki and /_/ = Ii, we need only 
consider I >0. Then from [27] one can derive the result, valid for Z = 0, 1, . . . 



/oo 
drr In rKf (mar) / Kf (ma) 



1 d rK 



2 dma V K 



[(-iy|l + 2E(-i)'-"fe]> 



2(ma) 



n=l 



K 



(5.42) 



where the modified Bessel functions on the right-hand side have argument ma. This gives 

-1 



d In dct 



Sch 



= — 2(ma)^ao 



ma 



dma \ Kq 



N 
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d 



ma 



dma 



In 



1 -1 



-2{may^ai 

1=1 



ma 



+ N 

d 
dma 



In 



H-N 



1 -1 



-N 



+ 



ma- m 



ma 



ma 



dma 

d 
dma 

d 
dma 



In 
In 



~1q 



-N 



+ N 



n -1 



N 



(5.43) 



The remainder of this section will be confined to an analysis of the asymptotic behaviour 
of detsc/i for large flux. The case of small fermion mass is also of interest as this limit will be 
controlled by the zero-energy bound states as noted in Sec. V.B. In addition, we conjecture 
that the low mass end of the integral in (2.1) will control the large flux growth of Indet^en- 
Hence we are led to consider the hmit N » 1 » ma of deigch- 



D. detsch for AT >> l >> ma 

The calculation of the above hmit requires the large I behaviour of ai in (5.42). Letting 
X — ma and using 



21 

X 



(5.44) 



we get 



1 d fKi_i\ I 



2 dx \ K] J X 



K, 



Ki 



i 



From [28] one finds for / >> 1 >> x 



1 1 o\ 1 / 1 13 



1 + - T^;" T + 



—x' + — + 



12 A J I V288 48 ?,2 J P \P 



(5.45) 



, (5.46) 



and hence 



Ki-i{x) ^ x_ 
Ki{x) ~ 21 



1 + - + 1 - -X 



1 o\ 1 



4 J I 



(5.47) 
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so that 



1 /I 



(5.48) 



Combining (5.48) with 



(5.49) 



vahd for Z, X > 0, we may conclude that 



lim Y^ai 



N-*oo 



1=1 



dx \ Ki 



0, 



(5.50) 



in which case the second, fifth and last terms in (5.43) vanish in the limit N — > oo. 



Considering the first term in (5.43) one finds that for a; << 1, 



ax \Ko 



In 



X 



-1 



+ 



In 



(f) 



-2 



(5.51) 



and 



X^Qq 



lull 



+ 



(5.52) 



so that the first term in (5.43) tends to — 1 for >> 1 >> x. 



Now consider the fourth term in (5.43). Combining [28] 



/ xe 



(5.53) 



for Z >> 1 >> X with (5.46) one finds 



X ^ \n[ "^'^^ 
dx \ Ki 



21 + N + O 



(5.54) 



which together with (5.48) gives 



l=N+l 



d , / h+N 
X— In 



dx \ Ki 



N 



-3 



=N+1 



(5.55) 
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For the range 1 < Z < A^, use (5.53) and K'jKi < for l,x,> to conclude, for >> 

i,;>i, 



dx V K] 



(5.56) 



Therefore, for x << 1 



oo 

lim a; 



— finite. 



(5.57) 



dx \ Ki 

Finally, consider the remaining third term in (5.43). For the range I > N the third 
term may be written 

-1 



N+l l=\ 

From (5.46) for N » 1 >> x,l>l 



d , Ii 



N 



l+N , 



(5.58) 



dx \Ki+nJ Ii \N + lJ \N + lJ' 



(5.59) 



since ///// > for /,a; > 0. Equation (5.59) together with (5.48) ensures that 
hm^v-^oo Eat+i = finite. 

For the range 1 < I < N let 

where we used /_/ = Ii. We find for x << 1 



bi = x^ai 



(5.60) 



bi = --lnx + 0(l) 



&2 = ^(l-^)+0([xlnxn 



The terms in (5.60) for 3 < Z < and x « \ behave as 



a, = i(/-2)--^ 



2/ -3 



8 {l-\f{l-2Y 



+ 0(x^). 



(5.61) 



(5.62) 



(5.63) 
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dx 



K, 



2{N - I + - 1) 

3NH -N^- AN^ - 31^ N + 8NI - SA^ 
(/ - iy{N - I + 1)2(/ - 2){N -1 + 2) 

+0(x^), 



(5.64) 



so that 



N 



^6, = -ln7V + 0(l) 

1=3 ^ 



(5.65) 



Combining the results obtained above on the sums in (5.43) together with (5.61) and (5.65) 
give 



d In detsc/i 



9$ 



In )+0(l). 



(5.66) 



The —\a.N growth is consistent with the data m Fig. 3. Recall that $ denotes — , the 

ZTT 



two-dimensional chiral anomaly, and that 



27r 



is the number of zero-energy bound states 



of the two-dimensional Pauli Hamiltonian (3.3). Assuming a smooth variation of det^c/i with 
we get, for ma « 1 << — 

ZTT 



\ndet sch = —(no. bound states)X 



, / no. l)()uii(l states\ , , 
In ( + 0(1) 



ma 



(5.67) 



We suspect that (5.67) holds more generally than for the finite flux magnetic fleld in (3.6). 



VI. ZERO MASS LIMIT OF detsch 

Consider either of the representations (2.5) or (4.19) for detsch- Does the m — limit 
of detsch exist? If it does, is it continuous in the sense that the massless Schwinger model's 
determinant is regained? We do not have any definite answers to these questions, but we 
suspect that they probably depend on whether the magnetic fields have zero or nonzero fiux 
as we will now explain. 
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Already at the level of second-order perturbation theory one runs into trouble for fields 
with $ 0. The determinant to order is 



Indetsc/. = -— / 77r4^|5(g)p / dz- 



27ri (27r)2 



z{l — z)q^ + m? 



(6.1) 



The finite flux field (3.6) gives 
In detsch = 



/e<l>\^ /-i^ ^ / a\m\ | [ 



m^O V27ry 
while the zero flux field 

B 



) ln(a|m|) + finite, 



(6.2) 



B(r) = - \5ir -a)- 5ir - 6)1 , 
r 



(6.3) 



gives 



Indetsc. = -{eBf t dzz{l - z) H dqq^^^ 
Jo Jo Z(l — 



[Joiqa) - Joiqb)]' 



z)q^ + m? ' 

which converges to the massless Schwinger model's determinant in the limit m = 0: 



(6.4) 



In deisch 



d\ \B{q)\^ 
27ri (27r)2 ^2 

= -(eB)^ \n(a/b), for a > & > 0. 



(6.5) 



Going beyond perturbation theory we have from (5.67) for the finite flux field (3.6) 

|e$ 



In det5c?i = 

m^O 27r 

|e$| >> I 



■ln(|m|a)-|- finite atm = 0. 



(6.6) 



(no. zero modes of H) x ln(|m|i?) + finite, 



(6.7) 



This result agrees with intuition, that is 

In = 

_det(p2 + ?77,2)J m — »■ 

where H is the Pauh Hamiltonian (3.3), and R is some natural length scale, such as the 
range of B. The massless Schwinger model's determinant is not regained. 



29 



Finally, Seller [12] defined the massless Schwinger model's determinant, det(l — e^^), by 
a renormalized determinant. This was done by making a formal similarity transformation 
and defining the determinant as det^en(l ~ ^K) with 

K{A) = (6.8) 

where K is considered as an operator on two-component square-integrable functions. The 
determinant det^en excludes the linearly divergent graph TrK but retains the graph TrK^, 
which is defined in some gauge invariant way. Assuming that the magnetic fields have zero 
flux and in particular that j (fr\A^\'^ < oo for all g > |, Seller was able to show that K E 
(the space of operators for which Tr{K^K)'^ < oo) for all n > 2 and that the spectrum of K 
consists only of the origin. This latter result is an expression of the triviahty of the massless 
Schwinger model and implies that all single-loop "photon-photon" scattering graphs of order 
and higher vanish [29]. These results allowed Seller to obtain the well-known result for 
the massless Schwinger model's determinant 

2 

lndet,en(l - eK) ^ J d''rA% (6.9) 
where ^4^ is in the Lorentz gauge. 

In the nonzero flux case falls off as ^, and K{A) ceases to be even a compact op- 
erator: the eigenvalues of K flU an open disc and K ^ Cn for any n > 1. Consequently, 
it is no longer possible to define the massless Schwinger model's determinant in terms of a 
renormalized determinant whose zeros refiect the eigenvalues of K{A). 

The lesson here is that the transition from a zero-flux magnetic field to one with nonzero 
fiux is not a smooth one and that the zero-mass limit of detsch will be fiux sensitive. 
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VII. NET RESULT 



Fermionic determinants are at the heart of fermionic field theories. In the case of QED 
the determinant in 3 + 1 dimensions in a static, unidirectional, finite (including zero) fiux 
magnetic field can be calculated from the determinant of the massive Euclidean Schwinger 
model, detsch, in the same magnetic field by integrating over the fermion's mass. Therefore, 
the massive Schwinger model is physically relevant. The calculation of det sch reduces to a 
problem in nonrelativistic, supersymmetric quantum mechanics, and the gauge invariance of 
det Sch is closely linked to the index theorem on a two-dimensional Euclidean manifold. The 
inclusion of mass quahtatively changes the determinant in 1 + 1 dimensions to the extent 
that the massless Schwinger model's contribution to detsch is cancelled by a contribution 
from the massive sector. Evidence was given that the zero-mass limit of det sch is not con- 
tinuous in the sense that the massless Schwinger model's determinant is not regained for 
nonzero fiux magnetic fields. 

It is believed that the first calculation of det5cfe for a finite fiux magnetic field is given 
in Sec.V. The behaviour of the determinant for large fiux and small mass suggests that 
the zero-energy bound states of the two-dimensional Pauli Hamiltonian are the controlling 
factor in the growth of In det 5cfc- If this is the case then the implication of this fact on the 
still unknown growth of the renormalized determinant of QED4 in the same magnetic field, 
which is determined by (2.1), remains to be seen. 
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FIGURE CAPTIONS 



FIG. 1. The logarithmically divergent contribution to the vacuum energy in QED2. 

FIG. 2. The contour of integration, C, in the /c-plane for the integral in (5.17). 

FIG. 3. Numerical calculation of the right-hand side of Eq. (5.41) for half-integral 
values of $ with ma = 1 and 10~^. 
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